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ON ' Using dimensional regularization, we compute the one-loop quantum and thermal 

■ corrections to the profile of the bosonic 1 + 1-dimensional 9s 4 kink, the sine-Gordon 

CN ■ kink and the CP 1 kink, and higher-dimensional ip 4 kink domain walls. Starting from 

the Heisenberg field equation in the presence of the nontrivial kink background we 
derive analytically results for the temperature-dependent mean field which display 
C^) ' the onset of the melting of kinks as the system is heated towards a symmetry restor- 

ing phase transition. The result is shown to simplify significantly when expressed 
in terms of a self-consistently defined thermal screening mass. In the case of do- 
main walls, we find infrared singularities in the kink profile, which corresponds to 
D I interface roughening depending on the system size. Finally we calculate the energy 

density profile of y? 4 kink domain walls and find that in contrast to the total surface 
tension the local distribution requires composite operator renormalization in 3+1 
. dimensions. 
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ON ! I. INTRODUCTION 

o 

Solitons [l| in scalar field theories in 1+1 dimensions, which we call generically kinks, 
have played an important role in increasing our understanding of various nontrivial aspects 
of quantum field theories, ranging from exactly solvable examples of strong/weak-coupling 
dualities 0, S] to the theory of topological defects generated at phase transitions with ap- 
plications in condensed matter physics [H as well as cosmology 

In a number of models one knows in closed form the spectrum of fluctuations about the 
kink background which allows one to perform complete calculations of one-loop corrections 
to the mass of the kinks. However, these calculations turn out to be full of subtleties, 
in particular (but not only) in the presence of fermions. For example, for the minimally 
supersymmetric kink a number of authors have concluded from explicit calculations that 
there was a cancellation of the one-loop effects on mass and central charge in a certain 
minimal renormalization scheme [fl, 0, 0], a result widely accepted since the mid 1980's. 
Only in 1997 two of the present authors have reopened this issue by demonstrating an 
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incompatibility of the methods employed for the supersymmetric kinks with known exact 
results for the nonsupersymmetric sine-Gordon model [9|, with correct results for the mass 
eventually being established in Refs. [ill . [l2| and a resulting puzzle concerning BPS 
saturation solved in Ref. 12j by the discovery of a new anomaly in the central charge 1 . More 



monopoles [14 
calculations 



14. 


15 




16 




17 



where a long-standing (20 years) but unnoticed discrepancy of direct 
with newer developments (notably Seiberg-Witten theory [l8j) was 
eventually cleared up. 2 

In the present paper we shall consider bosonic kinks at finite temperature 3 , applying 
and slightly generalizing the method used in Ref. [l2|, 21] to calculate the profile of kink 
to one-loop order. Following Ref. 2l| we start from the Heisenberg field equation. Other 
authors have used the method of the effective action with x-dependent background fields 
(see (IJ and references therein). Both methods are of course completely equivalent but we 
found the former to be simpler to implement. 

Quantum and thermal corrections to kinks have been considered in various approxima- 

24j . Here we shall 



tions of self-consistent, nonperturbative frameworks e.g. in Refs. [22|, [2c 
restrict ourselves to one-loop effects in a regime where perturbation theory is still reliable, 
refraining therefore from attempts to cover the physics of the symmetry restoring phase 
transition itself and the corresponding actual melting of kinks. However, we can reliably 
cover the onset of the melting of kinks and in this way provide benchmarks for more dar- 
ing approximations and approaches. Systematic calculations of one-loop corrections to kink 
profiles and domain walls at zero and finite temperature have been carried out before by 



several authors 25. l26l. 12 
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Our results are in agreement with those once the differ- 
ences in the renormalization schemes are taken into account. 4 However, we find significant 
simplifications in the corrections to the kink profile provided renormalization conditions are 
formulated in terms of a self-consistent thermal screening mass, which may be of practical 
importance in applications at non-infinitesimal coupling. 

We shall in turn cover three different models in 1+1 dimensions: the well-known exactly 
solvable sine-Gordon model with Z x Z 2 symmetry, a closely related massive CP 1 model with 
U(l) xZ 2 symmetry, and the familiar (p A model with Z 2 symmetry. All the above models have 
a discrete symmetry which is spontaneously broken at zero temperature and which leads to 
topological solitons (kinks). Since the y? 4 model is renormalizable in higher dimensions as 
well, we shall use it also to discuss domain walls in 2+1 and 3+1 dimensions. 

At sufficiently high temperature T one expects symmetry restoration and the disappear- 




1 Actually the conformal central chare 

2 For more extensive reviews see e.g. [19t |2C 

3 Finite temperature breaks supersymmetry so that supersymmetric kinks do no longer display features 
that are not also found in the bosonic case. Moreover, thermal contributions from fermions provide no 

iffjculty and are easily added on to what follows. 



special c 
4 In Ref. 



301 ] two of us have previously stated that we disagreed with Ref. 25| regarding results on the 



surface tension at zero temperature in a zero-momentum renormalization scheme, but agreed with Ref. 31 1 



As has been shown in Ref. 



29j, this apparent discrepancy was due to having compared two versions of 
The zero-momentum scheme considered in Refs 



30 



31] 



a zero-momentum renormalization scheme 
renormalizes the derivative of the two-point function with respect to momentum to one, while the scheme 
in Refs 



25 



291 ] does not introduce nontrivial wave-function renormalization, which is possible at one-loop 



order. Taking this difference into account resolves this issue. 
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ance ("melting") of the kinks. In 3+1 dimensions it is well known that perturbation theory 
allows one to derive the leading order result for the phase transition temperature T c , but 
in simple scalar field theories next-to-leading order results for the transition temperature as 
well as the order of the phase transition are beyond perturbation theory j32|, 33 . 



Therefore, let us explain the validity and limitations of perturbative methods in the 
various dimensions. In 3+1 dimensions, symmetry restoration in ip 4 theory with coupling 
constant A is brought about by a thermal mass term Am 2 ip 2 ~ XT 2 (p 2 which outweighs the 
wrong-sign mass term oc —fi 2 f 2 in the classical potential for sufficiently high temperature 
T > T c ~ [i/\/~X- After a resummation of the leading-order thermal mass, perturbation 
theory around the minimum of the effective potential has a loop expansion parameter ~ 
XT/m, where m is the (thermally corrected) mass of the fluctuations around the minimum, 
with m ~ \pl\x at low T in the broken phase, and m ~ Am for large T in the restored phase. 
With Ac 1, perturbation theory works for all temperatures except very close to T c where 
m gets parametrically small. (As long as m ~ /i, the expansion parameter XT/m < \f~X 
up to temperatures of the order of T c ; for T ^> T c , one has m ~ \f~\T and the expansion 
parameter is of order a/A throughout the symmetry-restored phase.) 

In lower dimensions, the situation is much more dire. In 1+1 dimensions, the coupling 
constant of y? 4 theory as well as of the sine-Gordon model has scaling dimension mass 
squared, and we have to assume that the loop expansion parameter at zero temperature 
X/m 2 <C 1. (For the CP 1 kink we have a dimensionless coupling constant replacing X/m 2 .) 
At finite temperature the expansion parameter is (X/m 2 ) x (T/m). High-temperature ther- 
mal mass terms are however only linear in T, Am 2 ~ XT/m. For symmetry restoration we 
would need |Am 2 | > m 2 , but this contradicts the requirement (X/m 2 ) x (T/m) <C 1. Hence, 
in perturbation theory we can reliably study the high-temperature limit T/m ^> 1 only as 
long as T/m <C m 2 /X, i.e. in the broken phase where |Am 2 | <C m 2 . It is therefore not 
mandatory to resum the thermal mass Am, although we shall find that it will be natural to 
do so. 

An important difference to the 3+1-dimensional case is, however, that Am 2 ~ XT/m 
in the 1+1-dimensional theory is generated only by Matsubara zero modes, whereas in the 
3+1-dimensional case the leading terms in Am 2 are generated by "hard", short- distance 
modes with wavelength ~ T" 1 , which for T/m 3> 1 are insensitive to the presence of a 
nontrivial kink background with inherent length scale ~ m~ x . 

In 2+1-dimensional theory, the situation is not better than in 1+1 dimensions. The loop 
expansion parameter at zero temperature is X/m <C 1. The thermal expansion parameter 
XT/m 2 <C 1 equally implies that the thermal mass squared Am 2 ~ AT <C m 2 , precluding 
a perturbative analysis of a symmetry-restoring phase transition. When going to high tem- 
peratures T/m ^> 1, perturbation theory is reliable only as long as T/m <C m/X, so again 
only the broken phase is accessible. 5 

In order to study the actual melting of kinks, nonperturbative methods are needed. (In 
the lower dimensional cases, also the symmetry-restored phase with T ^> T c is nonpertur- 
bative.) In some cases, other systematic expansions like large- N expansions j36l . 3?J can 



be put to work, but often (self-consistent) approximation schemes are employed which lack 
an expansion parameter controlling the approximations. In order to have credibility, such 



5 The critical temperature of the 2+1-dimcnsional (p 4 model has been estimated by renormalization-group 
methods in Ref. 34[ to be given by a relation of the form T c / [i oc fi/ Alog(cA/T c ) and the constant c 
subsequently measured on the lattice in Ref. 35 1. 
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approximations should be able to reproduce our perturbative results as limiting case. 

In Sect. |TT]we consider the sine-Gordon model and the massive CP 1 model, which have 
closely related fluctuation spectra, and we calculate the one-loop corrections to the field 
profile of the kinks in these models. In Sect. II I II we turn to the familiar p> 4 model, which 
has a more complicated fluctuation spectrum and correspondingly more complicated one- 
loop corrections. Since the latter model is renormalizable in higher dimensions, we shall 
use it to discuss one-loop corrections to domain walls in 2+1 and 3+1 dimensions, with a 
detailed discussion of our dimensional regularization and renormalization scheme, and its 
appropriate generalization at finite temperature. In particular, with the help of a remarkable 
identity for integrals involving the Bose-Einstein distribution function, we find that a self- 
consistent definition of the thermal screening mass removes certain artefacts in the one-loop 
kink profile. 

In the calculation of the profile of the higher-dimensional domain walls, we encounter 
infrared singularities associated with the massless modes that correspond to the transla- 
tional zero mode of the 1+1-dimensional kink. In accordance with Refs. (26|, 27, 29|, these 



singularities are interpreted as the field theoretic equivalent of system-size dependent in- 
terfacial roughening [38|, [39|, EL Finally we calculate the energy profile of the p> 4 kink 



and the corresponding domain walls and show that in contrast to the total mass (surface 
tension) the local energy density profile is ambiguous, depending on improvement terms to 
the stress tensor, and that in 3+1 dimensions composite operator renormalization through 
improvement terms is required. 



II. SINE-GORDON AND CP 1 KINKS 

We begin by calculating the one-loop quantum and thermal corrections to the field profile 
of the 1+1-dimensional sine-Gordon model and a massive version of the CP 1 model, which 
both happen to have a simpler fluctuation spectrum than the y? 4 model (in particular, no 
bound states). A full-fledged discussion of our method of dimensional regularization will be 
introduced in Sect. 1111} where we turn to the p> 4 model in 1+1 and higher dimensions, which 
also turns out to involve less trivial renormalization conditions. 



A. Sine-Gordon kink 



The Lagrangian of the sine-Gordon model is 

Tfl 4 

£ = --d^p + — 



COS I If 

m 



with a real scalar field <p. We have a discrete Z2 x Z symmetry, given by tp 
p — > tp + 2ir with bgZ. The field equation in 1+1 dimensions is 



m 



d^tp = + dl) p = -= sin — p . 



The constant solution of this equation yields the classical vacua 

m 



p = 2nn 



(i) 

-<p and 
(2) 

(3) 
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which break the discrete symmetry spontaneously. The kink solution interpolating between 
the vacua with n = (x = — oo) and n = 1 (x = +oo) is 



Am 

7t 



arctan e 



m(x— xo) 



(4) 



The energy of this configuration is M = 8m 3 /A. The equation for the fluctuations rj(x,t) 
around the kink solution, 



d^d^r](x, t) — m 2 r](x, t) cos 



Vx 



777 



-ip K {X) 



(5) 



becomes, using r)(x,t) = e lujt <f>(x) and the kink profile 



-d 2 + m 2 1 



(x) = OJ 2 (j)(x) 



(6) 



cosh rax / 

where we have set Xo = 0. From this equation we obtain the zero-mode with energy u = 0, 

'm 1 



0o (x) 



and the continuum with ujk = \Jk 2 + m 2 , 



2 coshmx 



777 



z tkx ( tanhmx - i— J 



(7) 



(8) 



We note the useful relation between the continuum and the zero-mode (Ref. (2l|, Eq. (12)) 

2m n , , 

(9) 



|^(x)| 2 = l-— 



We can now compute the correction to the kink profile in the presence of the fluctuations. 
To this end, following Ref. [2l[ , we interpret the equation of motion (j2j) as a Heisenberg field 
equation and write <p(x,t) — > ipk{ x ) + <p(x,t) = <-Pk{ x ) + 4>i{ x ) + v( x yt)- Here, 4>\{x) = 
{(p(x,t)), and the quantum fluctuation field r)(x,t) = f(x,t) — (pi(x) obeys (77) = 0. 

The equation of motion for 4>\{x) becomes 



o2 2 

c> — m cos 



Va 



77} 



-<Pk 



771 



sm 



771 



where 



(10) 



(11) 



(ri 2 )r C n = (V 2 }-(V 2 }\ x ^ 00 

is the renormalized propagator. To obtain Eq. (TlOl) from Eq. (j2]) we have employed the 
following renormalization. We have fixed 5m 2 and SX in the renormalization of the mass 
m 2 — > m\ = m 2 + 5m 2 and the coupling constant A — > X + 5X such that all one-loop (finite- 
temperature) graphs with one vertex cancel 41]. This condition implies that the ratio m 2 /X 
is unchanged under renormalization, and 



Sm 2 = m 2 y = ^ (V 2 )\,^ 



(12) 



6 



The propagator at finite temperature is given by 



ort = rE r = r gi±M-)|^, (13) 

where uo n = 2nirT are the bosonic Matsubara frequencies and n(u) = [exp(u/T) — 

is the Bose-Einstein distribution. Here we have dropped the zero mode of Eq. (171), which 

corresponds to the translational degree of freedom of the kink and which can be taken 



care of by collective quantization [42[]. This expression is UV divergent and thus requires 
regularization. In Sect. II I II we shall introduce our method of dimensional regularization 
adapted to solitons, but in this and the next section we shall suppress its details, since the 
finite results for the kink profile below do not depend on them. (If we were to calculate also 
local energy densities, as we shall do for the model of Sect. IIII[ we would need to be more 
careful.) The treatment of the zero mode in our way of dimensional regularization will also 
be discussed more fully in Sect. IHIt where we cover the more general case of domain walls. 

Simply subtracting Eq. (Tl3j) according to Eq. (TTIj) . and making use of Eq. (JSj) yields the 
renormalized propagator 



(V )ren = ~M X ) 



1 f°° dk n(uj k ) 

2m 



mn .Loo 2vr u\ 



(14) 



Inserting this result into the differential equation (|T0|) yields the correction to the kink profile 

, . . v^A ( \ 9 f°° dkn(oo k )\ sinhrax 

fa( x ) = — [— + m 2 ^ 5 . 15 

Yn ' 4m V2vr J^ 00 2tx w\ J cosh 2 mx V ; 

The total kink is then given by (Pk( x ) + 0i( a; ) with ^p K from Eq. (j4j). We see that there is a 
zero-temperature and a finite temperature correction, both being proportional to the deriva- 
tive of the zero mode. This is a consequence of the renormalization (fT2l) which subtracts the 
complete one-loop (seagull) diagram of the trivial sector including thermal contributions. 
Had we left out the latter (or any finite part), this would have produced extra contributions 
to 0i proportional to m-£^(p K = tp K + ^= 777,^/ cosh(ma;). 
With the large-temperature expansion 



2tt ojI 2m 4vr \T J ' v ; 



we see that the finite-temperature part of 4>i is suppressed compared to (fx by one power of 
the expansion parameter XT/m 3 (while the zero-temperature part is suppressed by A/m 2 ). 
Since this expansion parameter cannot be small in the case of symmetry restoration (see 
discussion in the introduction), this perturbative result is only valid for temperatures much 
smaller than the critical temperature of symmetry restoration. 

Identifying zero-temperature and thermal contributions of the mass counterterm 5m? = 
(^T=o m2 + drm 2 , we note that our renormalized mass m? = m\ — 5m 2 = (m 2 , — <5r=o m2 ) — 
Sj^m 2 = 77iy =0 + ?Ti j. differs from the renormalized mass at zero temperature by a negative 
thermal correction 

2 x 2 A f°° dkn(uj k ) XT 

m T = —o T m = / = h... tori ^> m. 17 

T 2 7_ 00 27r uj k Am K ' 
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FIG. 1: Finite-temperature corrections to the sine-Gordon kink (Pk( x ) + 4>i{ x ) according to the 
result (|15p for X/m 2 = 0.2 as a function of x times the zero-temperature mass itit=o and for three 
different temperatures, T/rriT=o = 0,10,20. The zero-temperature result is given by the dashed 
line, and with increasing temperature the kink becomes flatter. 



This means that thermal corrections tend to flatten out the potential by reducing the dif- 
ference between maxima and minima of the potential, which is proportional to m 2 (m 2 /A) 
with m 2 /X invariant. However, the distance between the minima and thus the value of (fx 
at x = ±oo remains fixed. In Fig. [T] we plot the kink profile for three different temperatures 
as a function of m^oi, showing that with increasing temperature the kink profile becomes 
flatter. (Had we plotted the profile as a function of mx with m the temperature-dependent 
mass, the kink would have appeared to become steeper instead.) 



The discussion of the finite-temperature kink in the mass-deformed CP 1 model is very 



similar to the sine-Gordon model of the previous subsection. The Lagrangian of the massive 
CP 1 model is 



with a dimensionless coupling r and a complex field 0. This model is renormalizable in 1+1 
dimensions and requires a coupling constant counterterm r — > r + 5r which is equivalent to 
wave function renormalization and a mass counterterm m 2 — > m 2 + 5m 2 . 

The Lagrangian (JTSJ) has a U(l) symmetry, — > e ia <f), and a discrete Z 2 symmetry, 
— > 1/0'. Since we work in 1+1 dimensions, the continuous symmetry cannot be sponta- 
neously broken j43|. A vacuum expectation value of the field rather breaks the discrete 
Z2 symmetry spontaneously. The classical potential 



B. CP 1 kink 




(18) 



V 




(19) 



is minimized at = and |0| = 00, i.e., by the south and north poles of the 2-sphere 
representing the compactified complex plane. (Both minima are invariant under U(l) as 
required) . 



s 



For a static solution, we can rewrite the Hamiltonian as 

r 



H = (d x <jt - m0t) (d x <f> - m0) - rmd x (l + 0*0) 1 (20) 



and read off the Bogomolnyi equation {d x — m)0 = 0. Thus the classical solution is 

tp K ( x ) = e ™(*-*o) ; (21) 

and the classical energy is M = rm. We again set Xq = 0. 

Next we determine the spectrum of the fluctuations in the presence of the kink. From 

t)C r 

g T = -[p(n-m 2 )0-20t((^) 2 -m 2 2 )], P=1 + 0V, (22) 

it is clear that the field equation of the fluctuations is obtained by setting = 77, keeping p 
fixed. This yields 

U - dl + ^k r md x + m^—^f) V {x, = 0. (23) 

With the ansatz 

V (x,t)=r- 1 /\l+tf < )g(x)e^ t (24) 

the field equation for g(x) becomes Eq. ([6]) of the sine-Gordon model. Consequently, as in 
the sine-Gordon model, we have a zero-mode 

, . [m 1 

^ ) = V^hw (25) 

and the continuous spectrum ujk = Vk 2 + m? with 

9k{x) = —e ikx ( tanhmx - i— ) . (26) 
u k \ mj 

To obtain the corrections to the kink profile induced by the fluctuations we write, analo- 
gous to above, 4>(x,t) = ipK(x)+(pi(x)+ri(x,t), and derive the equation for 0i (x). The terms 
with (77*77) are obtained by expanding the term — 2((p* / p)[2d x ipKd x r] — 2m 2 (pKf]} i n Eq. (|22|) 
to first order in 77*. It is obvious from (1221) that the propagator is diagonal: (rjr]) = 0. Hence 
we may set = <px + 4>i(x) + Vi( x i with real 0i(x) and 77! (x) and replace (77*77) by {rjifji). 
Taking into account also the mass counterterm 6 this yields 



1 + <Pk 1 + 1 Pk J 0- + ( PkF 1 + 1 Pk 

with 



-00 2vr 2^ 



6 The coupling constant counterterm does not contribute because it amounts to wave function renormal- 
ization and thus multiplies the classical field equation. 
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Using that (r]i(d x - m)^) = (\d x - m) (???) and {\d x - m) (l + (p 2 K ) 2 = m(<p 2 K -l)(l + <p 2 K ) 
we find that a mass counterterm of the form 



oo 



W = 4r"W/ ^ 1 + 2n( " fc) (29) 
2?r 2uj k y ' 



corresponds to a UV-subtracted quantity 

(vlU = r-\i + AT jT g 1 + ^ k) [\ 9k (x)\ 2 - 1] 



--^(i + ^(x)^- + 2m y-^j, (30) 

since the functions (^(x) and <7o(x) obey the same identity ([9]) as in the sine-Gordon case. 
With the explicit form of go(x), Eq. (125j) . we finally obtain a vanishing result for the right- 
hand side of the differential equation 



- rn)rji) vcn = - m^j (^i) rcn = , (31) 

leading to the remarkably simple result 

M*) = 0. (32) 

This means that the kink profile remains unaltered at finite temperature, provided the mass 
is renormalized according to fl29|) . The corresponding renormalization condition turns out to 
be an on-shell mass renormalization in the trivial sector including all thermal contributions 7 . 
The thermal correction to the mass has the same form as in the sine-Gordon model, with 
the replacement A — > 8r _1 m 2 in Eq. (|17[) . Since it is negative, this means that the kink 
e mx interpolating between the minima at <fi = and <fi = oo becomes spread out with 
increasing temperature. For the height of the potential between these minima, also the 
coupling constant renormalization is required. Direct calculation shows that the one-loop 
self-energy equals (p 2 + 3m 2 ) times a momentum independent expression, which implies that 
wave-function renormalization gives the counterterm 5r/r = ^5m 2 /m 2 . Hence, the thermal 
corrections of mass and wave function both work in the direction of diminishing the potential 
given in Eq. (fT9|) as the temperature is increased. 

Once again we can only determine the onset of the melting of the kink. The above 
discussion for the sine-Gordon model applies with r" 1 replacing the dimensionless coupling 
parameter X/m 2 . The requirement of a small expansion parameter at finite temperature, 
r _1 T/m 1, precludes the consideration of symmetry restoration as this would need 5m\ ^ 

2 



m 



7 In the CP 1 model the one-loop self-energy diagram is momentum-dependent, but it remains Lorentz- 
invariant at finite temperature, so that one does not have to distinguish between e.g. a screening mass or 
a plasmon mass. 
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III. ip A KINK AND DOMAIN WALLS 

The Lagrangian for the kink with \p 4 interaction is 

C = —d^p - ±(p 2 - v 2 ) 2 , t>WA, (33) 

with Z2 symmetry p — > —p. 

To one-loop order we introduce counterterms v 2 — > v 2 + 5v 2 and A — ► A + 5X and Sfi 2 = 
<5(At> 2 ) = A5f 2 +f 2 5A (leaving out wave-function renormalization which would be needed only 
at two- loop order in 3+1 dimensions). 5v 2 will be chosen such that in the topologically trivial 
sector tadpole contributions are subtracted completely (including thermal contributions). In 
the absence of wave-function renormalization, 5X is fixed by a renormalization condition for 
the mass of fluctuations in the trivial sector, which at tree level is m 2 = 2/i 2 . (Note however 
that in the spontaneously broken model (I33p the counterterms 5v 2 and 5X imply the mass 
counterterm 5m 2 = —\5v 2 + 2v 2 5\ ^ 25u 2 .) At zero temperature, the various possibilities 
have been discussed in detail in Ref. [30j. Below we shall address this question at finite 
temperature and single out one particularly natural renormalization condition. 

The equation of motion in the p> A model reads 

d^tp + pptp - Xp 3 = 0. (34) 
A kink at rest at x = x$ which interpolates between the two degenerate vacuum states 

V = ±-7= = ±v (35) 

is classically given by [l| 

Pk(x) = vt&nh (^{i(x — xq)/ V^j , (36) 

and its energy at tree-level is M = (v / 2/Lt) 3 /3A. From now on we set xq = without loss of 
generality. 

The kink can be trivially embedded in d+1 spacetime dimensions where it represents a 
domain wall separating the two vacua of the model. In this case, M has the meaning of 
a surface tension, i.e. energy per unit transverse volume. In the following we shall make 
d continuous and use this for dimensional regularization in our renormalization program 



301, 44, 45 



Fluctuations rj(x,t) about the classical kink solution are simple plane waves in the 
transverse directions y with transverse momentum £, i.e., r)(x,t) = |^e~ lw " f+rf ' 2 '0j(x), with 
u 2 e = uj 2 + £ 2 , where the x-dependent part <fri(x) is then given by the 1+1-dimensional 
fluctuation equation 

{-d 2 x - /i 2 + 3A0 2 ,) Ux) = ufyiix). (37) 

The spectrum of the 1+1-dimensional fluctuation equation [l[ has a zero-energy solution 
(w = 0), 

^W-^ir (38) 
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a bound state with energy ub = v / 3m/2, 



3m sinh^f 
~T cosh 2 ?f ' 



(39) 



and a continuous spectrum with energies ujk = a/A; 2 + m 2 , 



m 2 e lkx 



Auk-sjujl-io 



B L 



-3 tanh 2 h 1 + 4— + 6«— tanh 

2 m 2 m 2 



(40) 



Below we shall use the relation 21 



2m 



|0 fc (x)| 2 = l-^0 2 (x) 



■in 



(p 2 B (x) . 



(41) 



A. Kink profile 

Interpreting flMj) as a Heisenberg field equation and writing <p(x,t) = (Pk(x) + (pi( x ) + 
r)(x,t) with the quantum fluctuation field obeying {rj) = 0, we have the following equation 
for the one-loop correction to the kink profile, 



[d 2 x + fi 2 - 3Xip 2 K (x)] <j>i(x) = 3Xip K (x) 

= 3\<p K (x) (V 2 )ren.(x). (42) 



This equation is formally unchanged in the presence of a nonzero wave-function renormal- 
ization Z, although the counterterms Sv 2 and 5X will of course be different in different 
renormalization schemes. However, if two renormalization schemes, one with Z = 1 and one 
with Z 7^ 1, have the same renormalization condition for the mass m (e.g. that it should be 
the pole mass in the trivial sector without a kink), the one with nontrivial Z differs from 
the other one only by an extra contribution 8 A01 = — |(Z — 1 
Continuing with Z — 1, we have 

fa" = 3 <,%.., , =3 /^f 1+2 " K '' («) 

\ ' / I trivial sector J ( 2 7r) d 2uJ U V ' 

where we absorb the complete quantum and thermal contribution of the tadpole diagram 
in the renormalization of v (see Fig. [2^). Writing everything out for d+1 dimensions, this 



8 As one can check (cf. Sect. 2.2.2 of Ref. [30(), the net difference of the r.h.s. of ((42]) is then +X(Z - l)(fix, 
which leads to the extra contribution A<pi = ~\{Z — 1)<Pk- 
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leaves us with 



x 



d d - 1 £dkl + 2n(ou H ) 
(2n) d 2^ 



3m 2 /4 



1 



fc2 + cosh 2 (ma;/2) 
(3m 2 /4) 2 



d-i. 



(P + m 2 )(A; 2 + 
1 + 2n(ujBi) 3m / 



4 



-A/" d (T) 



(27r) d - 1 
3m 



2(jJb£ 

1 



<5A 



cosh 2 1 

2 

m 



mx/2) 
1 



cosh (mx/2) 



cosh (mx/2) 



(44) 



8 cosh 4 (mx/2) 6 A 2 cosh 2 (mx/2) ' 

where = V^ 2 + ^ 2 + ^ 2 and ujBt = \/^ 2 + 3m 2 /4. We have now also included the 
zero mode of the 1+1-dimensional fluctuation equation about the kink background, which 
in d > 1 spatial dimensions is a massless mode (the Goldstone mode associated with the 
spontaneous breaking of translational invariance), 



M d>x (T) 



d d ~H l + 2n(u M ) 



(45) 



where the T = part vanishes since in dimensional regularization scaleless integrals are 
identically zero. For d = 1 + e we discard the thermal part as well, since the latter is finite 
and we can set e = and uoe = there, assuming that this genuine zero mode is treated by 
collective quantization 42]. Thus we have 



Ad=i +e (T) = 0; N d (T = 0) = 0. 



(46) 



46 , and we 



In keeping a nontrivial thermal contribution J\fd for d > 2 we differ from Ref. 
shall later in Sect. IIII Bl provide evidence for the necessity of this. 

As mentioned above, in d = 3, 5X which appears in the last term of Eq. (1441) has to absorb 
a UV divergence of the momentum integrals, whereas in lower dimensions one could also be 
content with a minimal renormalization scheme, where only tadpoles are renormalized. In 
Refs. [H, S3 it was observed (for d = 1) that at zero temperature the requirement of an 
on-shell mass simplifies the result greatly - all terms proportional to 1/ cosh 2 (mx/2) then 
cancel. 

At nonzero temperature, there is actually not just one on-shell mass, but thermal masses 
are in general momentum dependent: at zero momentum, the thermal mass gives the plasma 
frequency above which propagating modes appear. Because of the absence of manifest 
Lorentz invariance, the effective mass of the propagating mode is generally not a constant. 
Thermal masses at frequencies below the plasma frequency give inverse spatial screening 
lengths. As we shall show now, it is the static screening mass which has to be employed in 
the renormalization condition in order that the above mentioned simplification occurs. 

Requiring that the renormalized parameter m be equal to the thermal static screening 
mass means that we have to subtract the self-energy diagram at zero frequency and imaginary 



spatial momentum q' 



-m 



The mass counterterm follows from (1331) after substituting 



9 In a nonabclian gauge theory, one finds that only this self-consistent definition leads to a gauge- 
independent result for the (Debye) screening mass 47] . In the present scalar theory, it will be seen 
to remove certain artefacts in the kink profile. 



13 



Y 



fiAz Sv 2 

+ X 



= 



(a) 



+ 



+ 



XSv 2 
--x 

-2v 2 5X 
--x 








(b) 



FIG. 2: Renormalization of the tadpole diagram (a) through the counterterm 5v 2 and of the self- 
energy diagrams (b) through the counterterm —5m 2 = X5v 2 — 2v 2 5X. Dashed lines correspond to 
the propagator in the trivial sector. 



<p — v + r), and reads 5m 2 = —X5v 2 + 2v 2 5X. As indicated in Fig. [2)3, the seagull diagram 
is already cancelled by the counterterm — X5v 2 , so the diagram with two propagators in 
the loop evaluated at q 2 = —m 2 defines 2v 2 5X. Using Feynman parametrization (which is 
straightforward to use in finite-temperature integrals at zero frequency, albeit not otherwise 
0) we find 



5X 



9X' 



dt 



d 



d-l. 



dk 



(27T) 



2n(u t ) 



where 



uj t = ^/k 2 + £ 2 + m 2 [l-t(l-t)]. 
In the T = part of 5X, one can easily integrate over the Feynman parameter to find 

6X ( T =o) _ 9X2 f d d ~ l ldk 1 



(27r) d ^k 2 + £ 2 + m 2 (k 2 + £ 2 + 3m 2 /4) ' 



(47) 
(48) 

(49) 



For later use we quote the closed-form result for this counterterm from Ref. [30 

,d-3 



a< T =°) = 9A 2 



m 



'3-(i^ 



(47T 



d+l 



(!) 



l± !V 2 ' 2' 2' 3'' 



(50) 



However, keeping the integral over transverse momenta I unevaluated and using 

1 2 



dk 



Vk 2 + a 2 (k 2 + b 2 ) bV~c 



b 2 



arctan 



b 2 



(51) 



one can verify that all zero-temperature terms proportional to 1/ cosh 2 (mx/2) in (|4*4|) cancel 
already under the integral over the transverse momentum I. Thus the cancellations observed 
in Refs. m, [m for zero temperature and d = 1 also take place for arbitrary d. 
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Turning to the thermal contributions one finds that the individual integrals involving the 
Bose-Einstein factor n cannot be evaluated in closed form. By numerical integrations one 
can however verify without difficulty the rather abstruse identity 



+ 



o 



1 , f dk 
dt / — 

2tt 



n(u t ) n'(ut) 



dkn(Vk 2 + £ 2 + 



2tt 



m 2 k 2 + ?f 



n( 



f 2 



\m 2 ) 



|m 2 



(52) 



where u>t was defined in Eq. (Hall . Using this to evaluate the Feynman parameter integral 
in SX, Eq. (14711 . and inserting the result into Eq. (14411 shows that at finite temperature it is 
the renormalization prescription of vanishing tadpoles and m being the finite-temperature 
screening mass which absorbs all terms proportional to 1/ cosh 2 (ma;/2) in (I44p . We then 
have 

1 



(r7 2 ) re n.(x) = -m d ~ l A d (T ', 'm) 



cosh (mx/2) 



(53) 



with 



A d (T/m) 



3m 



2-d 



d d ~H l + 2n(u Be ) 



3m 



2-d 



4 

g m 5-d 
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d d ~H 



2u>Be 
dk 1 + 2n(idke) 



-N d (T) 

1 



(27r) d 2uj m (k 2 + m 2 )(k 2 + ^f) 



(54) 



which when inserted into Eq. 



yields 



v m 



^A d (T/m) 



cosh (mx/2) 



v l m d ~ 



l A d {T/m) 



sinh(mx/2) 
cosh 3 (mx/2) 



(55) 



for the correction to the kink profile. In Fig. [3] the effect of such a correction to the classical 
kink profile is shown for (large) negative and positive coefficients. 

In a generic renormalization scheme, where (rj 2 ) rcn .( x ) h &s either uncancelled 
1/ cosh 2 (mx/2) terms or uncancelled constant terms (by incomplete tadpole subtraction), 
one would find additional terms oc mx/ cosh 2 (mx/2) in the kink profile correction 0i (aside 
from a different function A d (T/m)). A posteriori, the "bare" mx figuring in the latter term 
can be understood as an artefact of incomplete renormalization since the particular on-shell 
renormalization scheme considered above is evidently able to absorb these terms such that 
all mx appear only in exponentiated form. 

So the simple result (1531) depends on a renormalization scheme where m is the thermal 
static screening mass. This differs from the zero-temperature mass by the finite difference 



m — m 



T=0 



3A 



d 



d-it 



dk 



(27T) 



3m 2 /2 \ 9Xm f d d ~ x t n(uj B i) 



k 2 + 



(2vr) d - 1 u Bt 



(56) 



where we have again used the identity (152|) . 
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mx/2 

FIG. 3: The one-loop corrected ip 4 kink profile i^k{x) + </>i(x) oc tanh(z) +osinh(z)/ cosh 3 (2i) with 
z = mx/2 for a = —1 ... 1, showing the effect of positive vs. negative a oc up to nonperturba- 
tively large a. For d = 1, Ad is positive and increasing with temperature, as displayed in Fig. [H 
seemingly leading to a steeper kink at higher temperatures. However, this plot does not yet take 
into account that at increasing temperature the asymptotic values of the kink are reduced by the 
thermal part of 5v 2 , see Eq. (|43p . Including this effect, the kink profile does become flatter. 



1. 1+1 dimensional kink 

For d — 1, Md{T) is absent and all integrals in fl54l) are individually finite and the £- 
integration can in fact be dropped. The T = part is readily found to be 

M0) = ^ + l < 57) 



in agreement with Refs. |l2l. \2l\. The thermal part of A\{T) turns out to be strictly positive 
and growing linearly with T/m for large T/m. The full function is plotted in Fig. HI This 
corresponds to a positive parameter a in Fig. [3] which grows as the temperature is increased 
(with a having to remain sufficiently small so that perturbation theory is still valid). When 
plotted with fixed asymptotic values as in Fig. [31 the one-loop corrected kink profile which 
is slightly steeper than the classical one appears to become even steeper with increasing 
temperature. However, the asymptotic amplitude of the kink diminishes, because v 2 = 
Vq — 5v 2 = (vq — 5t=ov 2 ) — Stv 2 , and 5tv 2 , the thermal part of Eq. f|43l) . is positive and 
growing with temperature. Including this effect, the slope of the kink profile decreases at 
higher temperature. Moreover, also m has been renormalized such as to include thermal 
corrections and these reduce m as the phase transition is approached. When plotted in 
terms of a fixed zero-temperature mass rriT=o, the kink profile flattens even more quickly. 
At any rate, in perturbation theory we can only describe reliably the onset of the melting 
of the kink. 

The y9 4 kink in 1+1 dimensions has been studied extensively in nonperturbative self- 



consistent approximations to the two-particle irreducible effective action [22|, |23|, |24( . Com 



paring our Fig. [3] with Figs. 2 and 3 of [24j we find a somewhat different behavior even 



when the corrections are still small. This can be traced to the fact that in Refs. 22, 24 



variational ansatz for the dressed two-point function with coinciding arguments in the kink 
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FIG. 4: The function AdiT/m) in a doubly logarithmic plot for d = 1. This function determines 
the temperature dependence of the kink profile, see Eq. (|55l) . 



background, a nonperturbative generalization of our {f] 2 ) rcri ,{x), has been employed which 
is proportional to 1/ cosh 2 (mx/2). Our result f|53l) however suggests that a more adequate 
ansatz would involve 1/ cosh 4 (mx/2). In Ref. [22J the choice of the former was motivated by 
an analysis of the sine-Gordon model, where we have indeed obtained a renormalized two- 
point function proportional to 1/ cosh 2 (mx), see Eq. ( TT4~j) . But, as we have demonstrated, 
the generalization to the if 4 model is not justified. In this case, both the constant term in 
(rf) as well as the term proportional to 1/ cosh 2 (mx/2) are removable by renormalization, 
whereas in the sine-Gordon model only the constant term is. Note that in both cases the 
on-shell renormalized (rf) ren . to one-loop order turns out to be proportional to the zero-mode 
of that model squared. 



2. Domain wall profiles 



For d>2, the zero-temperature contributions in Eq. (l54j) are individually UV- divergent. 
Carrying out first the ^-integral in dimensional regularization we obtain (using formula 
(3.259.3) of Ref. gj}) 



4,(0) 



9 



1 



16 (47r) d / 2 

9 1 

16 (4vr) d / 2 



d 



2 

2- d 



d-4 

3 \ — 



1 

7T JO 

i r 



dx 



x 2 + \ 



TrT(^) 



l F 1 i 



d 1 6 
~'2 ;_ 



d 3, 
_; 4' 



(58) 



In the limit d — > 1 this of course reproduces Eq. (1571) . but for d — > 2 we encounter a 
singularity, because the square bracket in Eq. (158]) does not vanish for d — > 2. However, this 
divergence should not be an UV divergence, since, as we have shown above, mass, coupling 
constant as well as wave function renormalization can only modify constant terms and terms 
proportional to 1/ cosh 2 (mx/2) in (r] 2 ) and therefore not A^. Indeed, by taking into account 
A/2(0), which is zero in dimensional regularization, but whose integral representation is both 
UV and IR divergent, one finds that a small IR regulator mass £ 2 — > £ 2 + v 2 in (|45|) changes 
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(EHD to 



A d (0; v/m) 



9 



d 



16 (4vr) d / 2 V 2 
9 1 /2-d 



16 (47r) d / 2 



3 
4 

2 

3 ~ 3 



d-4 
2 



1/ 



d-2 
2 



> 2 + i)^ 

x 2 + \ 



0(d-2) 



which is finite for d — > 2 but now contains a term involving ln(z//m) 



A 2 (0) 



3 m 
m — , 

167T V 



(59) 



(60) 



so that the IR regulator cannot be removed. (The finite-T contributions in J\fd{T) have even 
more severe IR divergences.) Such IR divergences have in fact been discussed in the literature 
as being associated with the roughening of interfaces, leading to a logarithmic sensitivity 
of the width of interfacial profiles to the linear system size L. This phenomenon has been 



described in terms of the so-called capillary wave model [38|, |39|, |40j, with the capillary 
waves corresponding to the massless modes associated with the kink zero mode [26[. In a 
field theoretic treatment similar to ours, but with the interpretation of the Euclidean action 
of the (p A model as the Landau- Ginzburg Hamiltonian of statistical field theory, interface 
roughening has been discussed using cutoff regularization in Ref. (29|, and recently in d = 
3+e dimensions in Ref. [27j]. Our results for the one- loop profile, Eq. (1551) and (1601 . agree with 
Refs. [29| and 27], the latter upon identifying v ~ 1/L. (Ref. [27[ used a finite system with 
quadratic interface as IR regularization and also determined the resulting sublogarithmic 
contributions.) 

also involve IR-finite contributions of 
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The results for the profile given in Refs. [! 
the form mxj cosh 2 (ma;/2), which are due to the fact that there the results are expressed 
in terms of a static correlation length defined by a renormalization point q 2 = instead of 
q 2 = —m 2 . As we have shown, the latter definition, which is required to give m the meaning 
of an inverse exponential screening length, quite generally eliminates such contributions. 

For d = 3, i.e., in the 3+1 dimensional case, we find that the T = part is finite, yielding 
the (to our knowledge new) result 



A 3 (0) 



V3 



32vr 



(61) 



Now the thermal contributions are logarithmically IR divergent, because of J\fz(T), which 
leads to 

MT)~-4— (62) 
167r m v 

for an IR momentum cutoff v ~ 1/L, which indicates interfacial roughening effects also in 
the context of domain walls of relativistic field theories at finite temperature. 



B. Local energy density 



We now turn to the local energy density e(x), which, following Ref. [21J, we decompose 



as 



= ecas(£) + Ae Cas (x) + Ae (0l) (x). 



(63) 
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Here ecas( x ) represents the local energy density of a domain wall in d+1 dimensions (as- 
suming dimensional regularization) due to the sum over zero-point energies, whereas the 
contributions Aecas(x) and Ae^(x) are total derivatives which do not contribute to the 
integrated total energy (or surface tension) and which have been identified in Ref. [2l| . 

Subtracting off the energy density of the topologically trivial sector (including thermal 
contributions), the local energy density from the sum over zero-point energies is given by 



ecasUc) 



d 



d-l. 



(2 



7T 



-m 



id— 1 



[1 + 2n{u 0l )](t)l{x) 



d 



d-l. 



d d ~Hdk oj u 
{2n) d 2~ 



[1 + 2n(w w )] 



(27r) d - 1 2 
20 2 (^ 



[1 + 2n(u B£ )](f> 2 B (x) 



■~6v 2 [(f? K (x) 



5X r O 



■J' 



k 2 + m 2 



,21 2 



+ 



k 2 + ^f 



(64) 



where in the second line Eq. (HI!) has been used again to rewrite \<pk{x) \ 2 — 1 in terms of 4>o{x) 
and 4>b{x), and where we have included the contribution from the counterterms. (Recall 
that on he = Vk 2 + £ 2 + m 2 , ujbi= \J P 2 + 3m 2 / '4, and cj ^ = V^ 2 .) The first term in Eq. (j64"]) 
is the contribution from the massless modes corresponding to the zero mode of the 1 + 1- 
dimensional kink. Its zero-temperature part is eliminated by dimensional regularization, but 
its thermal contribution is that of black-body radiation in d — 1 spatial dimensions. It is 
to be omitted for d — 1, but contributes nontrivially in the case of domain walls. If it was 
not included (as done e.g. in Ref. [HI), the total one-loop surface tension J dxe^ would 
not vanish in the limit m — > 0, i.e. at the second-order phase transition where the kink has 
melted completely. To see this, note that for fixed i 3> k, one has IR singular limits m — ► 
of f dk( k2 2 m2 + fc 2 + ^2/ 4 ) in the second line of Eq. (164)) which cancel the explicit factor of m 
there. This yields thermal contributions which indeed compensate the m — *■ limits of the 
thermal contributions of the bound state modes and the massless modes. 

In contrast to the case of the kink profile we observe that in ec a s the massless modes 
do not lead to IR problems in either 2+1 or 3+1 dimensions, but we shall see that the IR 
singular kink profile plays a role in the remaining two (total derivative) contributions to 
e(x). The total energy (or surface tension) M is already determined by the above (IR-safe) 
expression, 

M= I dxe Cas {x). (65) 



(The integration over x is evaluated readily using J dx<p 



x) 



J dx(j>\ 



x 



V. 



For the local distribution of the energy, however, the total derivative terms Aecas(^) and 
Aeu 1 )(x) are relevant. The first of these comes from a surface term associated to a partial 



integration of the spatial gradients in the kinetic energy [21| and reads 



Cas 



4 x 



X 



5\ 



m 



6 A 2 cosh 2 (mx/2)_ 



(66) 



with (?7 2 ) r cn. defined in Eq. ( HSI) . 

The second contribution to the local energy density comes from the correction to the 
kink profile, </>i that we considered above. This also does not contribute to the total energy, 
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because the classical kink corresponds to a stationary point of the classical energy, but it 
gives a local modification of the energy density according to 



Ae 



(0i) 



d x {(f>id x tp K ) . 



(67) 



Inserting the results of Eq. ( |53|) and Eq. ( 155|) into these two expressions, we find that the 
terms with the function Ad, which appears in both, add with equal magnitude, yielding 



Ae Cas + Ae ( , 



= \m d+1 A d (T) 



~cosh 4 f£ 



+ 



cosh 6 ?f 



- SX^ [4>l{x) - 0|(x)] . (68) 



This localized contribution with zero total energy is finite in 1+1 dimensions, where it has 
been evaluated at zero temperature in Ref. [2lj], but in 2+1 and 3+1 dimensions it inherits 
the IR divergences found in the prefactor Ad of the one-loop kink profile discussed in the 
previous section. In 3+1 dimensions it is moreover UV divergent because of the appearance 
of 5 A in Ae Cas - 

Let us now check for UV divergences in Eq. 



), which was evaluated in |12|, [21| for 
the 1+1- dimensional kink, but only in integrated form for d+l- dimensional domain walls in 
Ref. 0- We shall now show that in the 2+1- dimensional case also the local energy density 
is rendered finite by on-shell renormalization (defined in the trivial sector). However, we 
shall find that in 3+1 dimensions there are uncancelled UV divergences in e c ~ (x) which 
require additional composite operator renormalization. 

Potential UV divergences can only come from the T = part (terms without n) in 
Eq. dSU). Using that cp 2 K {x) - v 2 = -fr(0| + 202) and [<f? K (x) - v 



212 



fp-0o we obtain 



-Cas 



X I 



X) - 

- ' 2 



+4>l{x) m 



d a 



(27T) d - 1 

d d ~ x ldk 



dk / 1 
oJbi + m I — — 

27T \ U) k i> 



(2tt) c 



1 



Integration over £ yields after some rearrangements 



k 2 + ?f 



m 



(69) 



C (T=0) 
-Cas 



X 



>B 



X 



(4vr) d / 2 2 



■(V)^ + m(l-d) 



— (k 2 + m 2 ) d ' 2 - 1 



+m 



[x 



^(k 2 + m 2 ) d / 2 -'J^4 iA 
2tt v k 2 + m 2 /4_ 

2\d/2-l 



m 



5A (T=0)^ 

6A 2 



(70) 



with 5A (T=0) given in fl50l) . This expression is finite for d — > 1, and includes an "anomalous" 
term proportional to (1 — d) J dk(k 2 +m 2 ) d ^ 2 ~ 1 , which would be missed in a naive momentum 
cutoff regularization. With it one recovers Eq. (18) of Ref. 2l|], except for the last term, 
which is due to 5X and which is finite by itself 10 . For d — > 2 the result is also finite (the 



10 The inclusion of this term corresponds to the finite renormalizations m — > m and A — > A in the final result 
(28) of Ref. [il | . 
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divergent T(— |) which is present in all terms except S A is cancelled after combining the UV 
finite expressions and using that J dk{k 2 + m 2 ) d l 2 ~ l is proportional to 1/T(1 — d/2)), but 
for d = 3 — 2e one finds divergent terms 1/e. All remaining finite integrals over k can be 
evaluated separately for d — 1,2,3 by using the substitution k/m = tant, but using the 
formulae given in Ref. [3(|, the following comparatively compact result can be derived, 



C (T=0), 
-Cas 1 



m 



2T 



3-d} 



X 



(47T 



(3/4)— /(d) 



(rf-i) 



+0 2 (x) d (3/4)" /(d) 



(71) 



where 



with the special cases 



f{d) = 2 Fi 



/(I) 
/(3 - 2e) 



3-d 1 3 



'2' 2' 



TT 



(72) 



2^3' 



/(2) 



^^3, 



TT 



V3 



-In 



(73) 



These can be used to write out closed-form results for the local Casimir energy density in 
1+1 dimensions (calculated previously in Ref. j2l|) and for the 2+1-dimensional domain 
wall. With J dx 0q = J dx <p 2 B = 1 one can readily verify that the total energy obtained 



from j dxe^ as °\x) agrees with the result given in Eq. (19) of Ref. 30|], which is finite for 



all d < 4. 

But for the local energy density we have UV divergent contributions for d = 3 — 2e, 



-Cas( X )|d=3-2e 



m 



6e (4tt) 



[0 2 (x) - <j> 2 B (x)] 



(74) 



which turns out to cancel incompletely with the UV divergence found above in Eq. (1681) . 



x)\ 



m 



d=3-2e 



2e (4?r) 



x) 



X 



(75) 



However, it is well known that the local energy density is ambiguous up to improvement 
terms to the energy-momentum tensor, and that an unim pro ved ener gy-m omentum tensor 



in general needs composite operator renormalization fed,|5i|, [52|, |53|, |54j]. This ambiguity 
signals that the local energy density is not a well-defined physical observable as such. If 
one defines the energy density as the 00-component of the gravitational energy-momentum 
tensor (variation of the Lagrangian with respect to the metric) then it will depend on the 
particular coupling to the gravitational field, which can also include a term in the Lagrangian 
of the form — 1£ y/—g(x)R(x)ip(x) 2 where R is the Riemann scalar. In a flat spacetime, this 
produces an additional contribution to the local energy density given by 



(76) 
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TABLE I: Divergences left after standard renormalization in the field profile </>i(x), the local energy 
density e(x), and total energy (mass or surface tension) for the various dimensions. "IR" denotes 
IR divergences which introduce a dependence on system size; "UV" denotes divergences requiring 
composite operator renormalization through improvement terms for the energy momentum tensor. 



dimension 


profile 4>\{x) 


energy density e(x) 


total energy M 


1+1 


finite 


finite 


finite 


2+1 


IR (all T) 


IR (all T, £ + \) 


finite 


3+1 


IR (T > 0) 


IR (T > 0, £ + ±), UV (£ + \) 


finite 



which can contribute already at tree- level for a "nonminimal" coupling £ ^ 0. 11 Even when 
£ = at tree level, as we have implicitly been assuming so far, quantum contributions will 
in general produce such a term. Indeed, the combination 0q — <fi 2 B appearing in the divergent 
contributions Eqs. (174)) and (175|) is proportional to d 2 ip 2 K (x), 

<Pl-<Pl = ^d^ K (x) (77) 

and thus can be removed by a counterterm —\5£ \f— g(x)R(x)(p(x) 2 with divergent 5£. 

Alternatively, having a nonzero £ already at tree-level produces further divergences of the 
form of (!74|) and (!75|) . namely 

Ae ? = -^ 2 )(x) = -4£Ae Cas . (78) 

This is UV divergent in 3+1 dimensions and also generally IR divergent because it involves 
the kink profile. By choosing £ = 1/2 we can in fact make sure that the correspondingly 
modified local energy density is independent of Ad(T) and thus always IR finite (but not 
UV finite in 3+1 dimensions). Alternatively, £ = 1/6 leads to a UV finite (but generally IR 
divergent) one-loop local energy density for d = 3 without infinite renormalization of £. (In 
a massless theory, £ = preserves conformal invariance in a gravitational background and 
thus £ = +1/6 is the value which corresponds to an improved stress tensor that is finite in 
3+1 dimensions.) 



IV. CONCLUSION 

We have calculated one-loop corrections to the profile of sine-Gordon and CP 1 kinks 
and y9 4 domain walls at zero and finite temperature. Using dimensional regularization, we 
have reproduced results for the one-loop field profile of 1+1-dimensional (bosonic) kinks of 



Ref . [12j, l21| , and have extended them to include thermal contributions, and also to cover 
higher- dimensional ip 4 kink domain walls. We have shown that a renormalization condition 
which defines thermal screening masses self-consistently at negative wave vector squared 



11 The associated integrated total energy vanishes for the kink, but in solitons involving long-range fields 
such as magnetic monopoles improvement terms to the energy-momentum tensor can also contribute to 



the total energy, see Ref. 55| . 
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simplifies the result and removes certain artefacts in the resulting one-loop corrected kink 
profile. 

However, in the case of domain walls, we have encountered divergences in local quantities 
that are not taken care of by standard renormalization of the parameters in the Lagrangian 
as summarized in Table [H On the one hand, we have found infrared singularities caused 
by the massless modes which are the higher- dimensional analogs of the translational zero 
mode of the 1+1-dimensional kink. In the 2+1- dimensional domain wall, these infrared 
singularities lead to a logarithmic sensitivity of the correction to the classical kink profile on 
the system size. This corresponds to the phenomenon of interfacial roughening in statistical 
physics, where the same model has been studied in a 3-dimensional Euclidean setting [27| . In 
3+1 dimensions, the one- loop kink profile at zero temperature turns out to be infrared finite 
in dimensional regularization, whereas thermal contributions lead to logarithmic infrared 
singularities, thus exhibiting the phenomenon of interface roughening in the context of a 
relativistic field theory at finite temperature. 

In the case of the (p 4 kink and the corresponding domain walls, we have also calculated 
the local energy profile and discussed its ambiguities. Depending on the underlying energy- 
momentum tensor, in particular the parameter £ in a possible —\^^—gRip 2 term in the 
Lagrangian, we generally found both infrared and ultraviolet divergences. In 2+1 and 3+1 
dimensions, the infrared divergences are the same as those found in the field profile, although 
for the choice C — o they can be eliminated from the energy profile. Ultraviolet divergences 
arise in 3+1 dimensions whenever £ ^ |, corresponding to an unimproved energy-momentum 
tensor, and these divergences can then be cancelled by improvement terms. However, both 
infrared and ultraviolet divergences drop out in the integrated energy density. Considering 
the manifestly finite thermal corrections to the energy density we have found that the con- 
tributions from the massless modes are crucial in ensuring that the energy density vanishes 
in the limit m — > 0, i.e. when the domain wall has melted in a second-order phase transition. 
In the case of the field profile, excluding these massless modes (as done e.g. in Ref. [46^ 1 
would have resulted in non-renormalizable ultraviolet-divergences in place of the (physical) 
infrared divergences. 
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